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ABSTRACT 

We present the results on the study of the global and local Wt-Z relation based on the first data available from the CALIFA survey (150 galaxies). 

This survey provides integral field spectroscopy of the complete optical extent of each galaxy (up to 2-3 effective radii), with enough resolution to 

separate individual Hn regions and/or aggregations. Nearly ~3000 individual Hn regions have been detected. The spectra cover the wavelength 
C* . ' range between [OII]3727 and [SIIJ6731, with a sufficient signal-to-noise to derive the oxygen abundance and star-formation rate associated with 
Q ' each region. In addition, we have computed the integrated and spatially resolved stellar masses (and surface densities), based on SDSS photometric 
I data. We explore the relations between the stellar mass, oxygen abundance and star-formation rate using this dataset. 
O . We derive a tight relation between the integrated stellar mass and the gas-phase abundance, with a dispersion smaller than the one already reported 
l *j ' in the literature (a" Alog ( /H) =0.07 dex). Indeed, this dispersion is only slightly larger than the typical error derived for our oxygen abundances. 
CZ3 ' However, we do not find any secondary relation with the star-formation rate, other than the one induced due to the primary relation of this quantity 
Cu | with the stellar mass. The analysis for our sample of ~3000 individual H n regions confirm (i) the existence of a local mass-metallicity relation 

and (ii) the lack of a secondary relation with the star-formation rate. The same analysis is done for the specific star-formation rate, with similar 
! results. 

Our results agree with the scenario in which gas recycling in galaxies, both locally and globally, is much faster than other typical timescales, 

like that of gas accretion by inflow and/or metal loss due to outflows. In essence, late-type/disk dominated galaxies seem to be in a quasi-steady 

situation, with a behavior similar to the one expected from an instantaneous recycling/closed-box model. 
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C) 1. Introduction mass-metallicity (A4-Z) relation is consistent with more mas- 
sive galaxies bei ng more metal-enriche d. It was confirmed ob- 

£> ' Metals form in stars as a by-product of the thermonuclear reac- servationally by iTremonti et all (I20041 hereafter T04), who 

tions that are the central engine of stellar activity. Once they have found a tight correlation spanning over 3 orders of magnitude 

completed their life cycle, stars eject metals into the interstellar in mass and a factor of 10 in metallicity, using a large sample 

medium, polluting the gas, which is the fuel for the new genera- of star-forming galaxies up to z ~ 0. 1 from the Sloan Digital 

tion of stars. Therefore, the star-formation rate (SFR), the stellar Sky Survey (SPSS). The At-Z relation app ears to be indepen- 

mass, the metal content and the overall star- formation history of dent of large-scale ( Mouhci ne et al.l [2007) and local environ- 

galaxies are strongly interconnected quantities. A fundamental ment (Hughes et al. ;2012l) . although it has been observed that 

open question in our understanding of galaxy evolution is based in high densit y environments metallicities are h igher than ex- 

on the details of these interconnections (are they local or global? pected (e.g. Mateus et al. 20071: iPetropoulou et alJl2012h . and in 

are both intristic or a consequence of the evolution?), and their lower density environment s satellities have higher me tallicities 

dependence on other properties of galaxies (are they affected by than central galaxies (e.g. iPasquali et al.l 120 1 OJ. |2012|) . Finally, 

the environment and how? how does the merging history of the the rel ation has been confirmed at a ll accessible redshift s (e.g. 

small interactions affect them?). ISavaglio etal . 2005; Erb et al. 2006; Maiolino lit~al"1l2008l) . 

The existence of a strong correlation between stellar Considerable work has been devoted to understanding the 

mass and gas-phase metallicity in galaxies is well known physical mechanisms underlying the At-Z relation. The pro- 

dLequeux et al.l Tl97 9: Skillman 1992). These parameters are posed scenarios to explain its origin can be broadly cate- 

two of the most fundamental physical properties of galaxies, gorized as: 1) the l oss of enriched gas by outflows (T04; 

both directly related to the process of galaxy evolution. The IKobavashi et alJ2007l) : 2) the accretion of pristine gas by inflows 
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(Finlator & Dave 2008); 3) variations of the initial mass function 
with galaxy mass (Koppen et al. 2007); 4) selective star forma- 
tion efficiency or downs i zing (Brooks et al. 2007; Elliso n et al.l 
120081: ICalura et al.lf2009l: IVale Asari et al.ll2009h : or a combina- 
tion of them. 

Despite the local nature of the star-formation processes, it 
has only recently become possible to analyze the M-Z rela- 
tion for spatially resolved, exte rnal galaxies. Early prece dents 
with more limited data i nclud e [Edmunds & Pagel ( 1984J) and 
IVila-Costas & Edmunds; dl992l) . who noticed a correlation be- 
tween m ass surface density and gas metallicity in a number of 
galaxies. Mor an et al.l (1201 2|) recently reported a correlation be- 
tween the local stellar mass density and the metallicity valid 
across all the galaxies in their sample. This result was derived 
analyzing the individual Hn regions of a sample of 174 star- 
forming galaxies, base d on slit spectroscopy. We in dependently 
confirmed this relation (Rosales-Orteg a et al.ll2.012h . using a sta- 
tistically complete sample of ~2500 H n regions extracted from 
a sample of 38 galaxies, us ing integral field spectroscopy (IFS), 
described in Sanchez et al.l (12012b). 

This local relation is a scaled version of the M-Z rela- 
tion. This new relation is explained as a simple effect of the 
inside-out mass/metallicity growt h that dominates the secular 
evolution of late-type galax ies (e.g. lMatteucci & Francoislll989l: 
Boissier & Prantzoslll999l) . combined with the fact that more 
massive regions form stars faster (i.e. at higher SFRs), thus 
earlier in cosmological times, which can be considered a lo- 
cal downsizing effect, simila r to the one observed in individual 
galaxies (e.g. Perez-Gonzalez et al. 2008). This explanation does 
not require a strong effect of inflows/outflows in shaping the M- 
Z relation, which can be naturally explained by secular evolution 
processes. 

Recently, evidence has been reported for a dependence of 
the M-Z relation on the SFR. Different authors found that these 
three quantities define either a surface or a plane in the cor- 
responding 3D space, i.e., the so-called Fundamental Mass- 
Metallicity relation (FM R) and/or Al-Z-Fundamental Plane, de- 
pending on the auth ors (lLara-Lopez et al.ll2010l: IMannucci et alJ 
2010; Yates et al. 2012). R egardl ess of the actual shape of this 
relation (lLara-Lopez et al.l 1201 2l) . its existence is still contro- 
versial. The described relationship implies that (i) for the same 
mass, galaxies with stronger star-formation rates, have lower 
metallicities, and (ii) for low-mass galaxies, the dependence on 
the SFR is stronger. Metallicity is a parameter that depends on (i) 
the star-formation (which enriches galaxies), (ii) inflows (which 
dilute and reduce the metallicity), and (iii) outflows (which eject 
metals out of the galaxy). Thus, the reported relation imposes 
restrictions to the ratio between the chemical enrichment, pro- 
duced by the SFR, and dyna mical timescales which reg ulates the 
dilution due to inflows (e.g. Ouillen & Bland-Hawthorn 2008). 
On the other hand, it also restricts the dependence between the 
amou nt of metals ejected b y an outflow and the strength of the 
SFR (Mannuc ci et al.ll2010l and references therein). In a stable 
situation, described by a simple instantaneous recycling model, 
no dependence on the SFR is expected. 

So far, all the studies reporting a dependence between the M- 
Z relation and the SFR are based on single aperture spe ctroscopic 
data, like the one provided by the SDSS survey (lYork et al.l 
2000). In spite of the impressive dataset provided by these sur- 
veys, with tens or hundreds of thousands of individual measure- 
ments, they present some well-known drawbacks: (i) single aper- 
ture spectroscopic data have a fundamental aperture bias, which 
restricts the derived parameters to different scale-len gths of the 
galaxies at different redshifts (e.g. lEllis et al.l 12005). This is a 



strong limitation for these studies, considering that both the oxy- 
gen abundance and the star-formation rate shows strong gradi- 
ents within galaxies (e.g. Sanchez et al. 2012b). Therefore, the 
derived parameters may not be representative of the integrated 
(or characteristic) ones; (ii) The fact that most of these surveys 
sample galaxies at a wide range of redshifts implies that the de- 
rived parameters correspond to different physical scales, for sim- 
ilar galaxies at different redshift, due to an aperture bias. This 
may induce secondary correlations difficult to address, in partic- 
ular, if some of these parameters present a cosmological evolu- 
tion and/or spatial gradients within each galaxy; (iii) while the 
spectroscopic information (oxygen abundance and SFR) is de- 
rived from this aperture limited dataset, the third analyzed pa- 
rameter, the mass, is derived mostly based on integrated photo- 
metric data. Moreover, in some cases it is derived using the infor- 
mation comprised in the spec troscopic data to estimate the cor- 
responding M/L-ratio (Man nucci et al.l 120101; iKauffman n et al.l 
2003a). Contrary to what is often claimed, we recall that these 
systematic effects cannot be reduced by the large statistical num- 
ber of data encompassed in the considered dataset. 

Some of these issues are partially solved when a large 
aperture is used to derive the spectra of each galaxy, as in 
the case of drift-sc an observations (e.g. iMoustakas et al.| [2010; 
Hugh es et al.l 120 12!). In this case, the derived abundance (basi- 
cally luminosity weighted across the optical extent of the galax- 
ies) has a better correspondence with the characteristic one (i.e., 
the value at 0.4p25, whe re pjs is the radius at a surfa ce bright- 
ness of 25 mag/a rcsec 2 . IZaritskv et al.ll 19941: lGarnettll2002l) . as 
demonstrated by Moustakas & Kennicutt (2006). On the other 
hand, the SFR can be directly derived from the integrated Ha 
emission across the considered aperture. However, these obser- 
vations have also some limitations: (i) in many cases they do 
not cover the entire optical extent of the galaxies, and the cov- 
ered fraction is different from galaxy to gala xy, which introd uces 
again an aperture uncertainty (Fig. 3 of San chez et al.ll201 U for 
an example); (ii) due to the lack of spatially resolved spectro- 
scopic information, they mix regions with different ionization 
properties and/or ionization sources, that could change substan- 
tially jrtdifferentjocations even within quiescent spiral galaxies 
(e.g. lSanchez et al] |2012b). This introduces uncertainties in both 
the derived abundances and the estima ted SFR . 

In spite of these limitations, Hughe set al.l (1201 2l) explored in 
a recent study the dependence of the Al-Z relation on the SFR, 
using drift-scan observations of 135 nearby late-type galaxies 
in the Local Universe. As expected, their dispersion around the 
M-Z relation is similar or even slightly higher than the one re- 
ported by T04. Although, th ey cannot reproduce th e results by 
IMannucci et al.l (12010) and Lara-Lope z et alJ (1201 Ol) . and in fact 
the scatter increases when they try to introduce a dependence on 
the SFR, they found a strong depedence of the metallicity with 
the gas fraction. This dependece could induce a correlation with 
the SFR in a large samp le of galaxies. 

M ore recently, iPerez-Montero et all (1201 2l) and ICresci et aD 
d2012l) studied the M-Z relation in a wide range of redshifts 
both using the zCOSMOS data, and the SDSS ones in the first 
case. They found either a negative trend of the metallicity with 
the SFR for a fixed galaxy mass, in the first case, or a dis- 
tribution of masses, metallicies and SFRs con sisten with the 
FMR, in th e lat ter one, both in agreement with IMannucci et al.l 
(2010) and Lara-Lopez et al. (2010). Despite of this agreement, 
IPerez-Montero et all (120121) noted that the dispersion around the 
M-Z relation is only reduced by 0.01 dex at all stellar masses 
when a secondary relation with the SFR is introduced, for an 
initial dispersion of ~0. 1 dex. As a comparison we should note 
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Fig. 1. To/? panels: IFS-based Ha maps, in units of 10 erg s _1 cm -2 arcsec" 2 , derived for two representative galaxies of the 
sample (color images), together with the detected H n regions shown as black segmented contours. Bottom panels: Radial distri- 
bution of the oxygen abundance derived for the individual Hn regions with abundace errors below 0.15 dex, as a function of the 
deprojected distance (i.e., corrected for inclination), normalized to the effective radius, for the same galaxies. The size of the circles 
is proportional to the He intensity. 



that th e decrease of the dispersion reported by Mannucci et al. 
(2010) was nearly half of th e original value, when introducing 
the same relation. Perez -Montero et alj (1201 2b considered that 
the lack of extreme SFR values in their SDSS galaxies may ex- 
plain this effect. 

At high e r red shifts the situation is more complicated. 
IWuvts et all d2012l) found that the distribution of their data 
is in agreement with the local FMR relation presented by 
iMannucci et alJ (l2010h . However, they do not see a correlation 
between metallicity and SFR at a fixed mass. On the contrary, 
Richard et al. (2011) interpret the largest metallicities found at 
the low mass range as a consequence of the FMR in combina- 
tion with a lower star - formation rate of these objects. In a similar 
way, iNakaiima et al.1 (1201 2l) show that the lower limit to the av- 
erage metallicities found for a sample of Lya emitters at z ~ 2.2, 
after stacking the individual spectra, is roughly consistent with 
the FMR relation. Finally, Magrini et al.l d2012l) . support the idea 
that the scaling relation between the mass, metallicity and SFR 
has a different origin (and shape) for "active" starbursts, more 
common at high redshifts, and quiescent galaxies, more frequent 
at lower ones. The fact that in most cases only the N2 indicator is 
accesible at high-redshift may introduce biases difficult to quan- 
tify. 



In order to explore the A4-Z relation in the Local Universe, 
and to bring light to its dependence on the SFR, while mini- 
mizing the bias effects previously described, we have used the 
integral -field spectroscopic (I FS) data provided by the CALIFA 
survey (Sanch ez et al.ll20 12alH CALIFA is an on-going explo- 
ration of the spatially resolved spectroscopic properties of galax- 
ies in the Local Universe (z <0.03) using wide-field IFS to 
cover the full optical extent (up to ~2.5 r e ) of ~600 galaxies 
of any morphological type, distributed across the entire Color- 
Magnitude diagram (Walcher et al., in prep.), and sampling the 
wavelength range 3650-7500A. So far, the survey has completed 
~l/3 of its observations, and the first data r elease, compris- 
ing 1 00 galaxies has been delivered recently (Husemann et al.l 
l2013h . 

The layout of this article is as follows: in Sec. [2] we sum- 
marize the main properties of the sample and data used in this 
study; in Sec. [3] we present the main analysis, and the derivation 
of the parameters analyzed in the article: the mass, metallicity 
and SFR for each individual galaxy and H n region; the global 
M-Z relation and its dependence with the SFR is explored in Sec. 
I3.2l and [3~3l a similar analysis for the local M-Z (or E-Z) relation 
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Fig. 2. Distribution of the stellar masses derived using the average M/L ratios described in Bell & de Jong (2001 ) and the B and V- 
band photometry (Y axis in all the panels) along the stellar masses derived using other procedures (X axis in all the panels): Top-Left 
panels Stellar masses derived using all the five-bands SDSS photometry and the M/L ratio provided by the SED analysis performed 
using PARAD ISE; Top-Right panel: Stellar masses derived using the 2MASS K-band photometry and the M/L ratio provided by 
lLonghetti & Sa racco (2009); Bottom-left panel: Stellar masses derived using the CALIFA datacubes spectrophotometry and the M/L 
ratio provided by PyCASSO/Starlight; and Bottom-right panel: Stellar Masses derived from CALIFA datacubes spectrophotometry 
and the M/L ratio provided by FIT3D 



is presented in Sec. 13.41 we explore the possible dependence on 
the specific star-formation rate (sSFR) in Sec. 13.51 our results 
and conclusions are summarized in Sec.|4] 



2. Sample of galaxies and data 

The galaxies under study have been selected from the CALIFA 
observed sample 1 . Since CALIFA is an ongoing survey, whose 
observations are scheduled on a monthly basis (i.e., dark nights), 
the list of objects increases regularly. The current results are 
based on the 150 galaxies observed using the low -resolution 
setup until July 2012. At that point, most of the color-magnitude 
diagram had been sampled by the survey with at least one or 
two targets per bin of magnitude and color, including galaxies of 
any morphological type. The CALIFA mother sample becomes 
incomplete below M r > -19 mag, which corresponds to a stel- 
lar mass of ~10 9 ' 5 M Q . Therefore, it does not sample low-mass 
and/or dwarf galaxies. Most of these galaxies are part of the 1st 
CALIFA Data Release lHusemann et alJ (l2013h . and therefore the 



datacube are accesible from the DR1 webpagqj Table [TJ shows 
the list of the galaxies analyzed in the current study, including, 
for each galaxy: (i) its name, (ii) redshift, (iii) V-band magni- 
tudes and B - V color and visual morphological classification 
(extracted from Walcher et al. in prep). In addition, we include 
the main properties derived along this article, as explained later: 
(i) the integrated stellar mass, (ii) characteristic oxygen abun- 
dance and (iii) the integrated star-formation. 

The details of the survey, sample, observation al strategy and 
reduction are explained in Sanchez et al. (2012a). All galax- 
ies were obser ved using P MAS (Roth et al. 2005) in the PPAK 
configuration dKelz et alJl2006h . covering a hexagonal field-of- 
view of 74"x64", sufficient to cover the full optical extent of 
the galaxies up to 2-3 effective radii. This is possible due to 
the diameter selection of the sample (Walcher et al., in prep.). 
The observing strategy guarantees a complete coverage of the 
FoV, with a final spatial resolution of FWHM~3", i.e., ~1 kpc 
at the average redshift of the survey. The sampled wavelength 
range and spectroscopic resolution (3745-7500A, A/AA ~850, 
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Fig. 3. Distribution of the oxygen abundances estimated for each 
galaxy at the effective radius (Re), using the radial gradients, and 
the oxygen abundances derived directly from the spectra inte- 
grated over the entire field-of-view (int). 



for the low-resolution setup) are more than sufficient to ex- 
plore the most prominent ionized gas emission lines, from 
[OII]/B727 to [SII]^6731, on one hand, and to deblend and sub- 
tract the underlying s tellar population, o n the other hand (e.g. 
ISanchez et aDl2012at iKehrig et alJl2012l Cid Fernandes et al., 
submitted). The dataset was reduced using version 1.3c of the 
CALIFA pipeline, whose modifications with respect to the one 
presented in [S anchez et al. (2012a) are described in detail in 
Huseman n et alj d2013l) . In summary, the data fulfill the fore- 
seen quality control requirements, with a spectrophotometric 
accuracy better than a 15% everywhere within the wavelength 
range, both absolute and relative, and a depth that allows us 
to detect emission lines in individual Hn regions as weak as 
-lO^ergs- 1 cm" 2 A" 1 , with a S/N-3-5. 

The final product of the data-reduction is a regular-grid dat- 
acube, with x and y coordinates indicating the right-ascension 
and declination of the target and z being a common step in wave- 
length. The CALIFA pipeline also provides the propagated error 
cube, a proper mask cube of bad pixels, and a prescription of how 
to handle the errors when doing spatial binning (due to covari- 
ance between adjacent pixels after image reconstruction). These 
datacubes, together with the ancillary data described in Walcher 
et al. (in preparation), are the basic starting point of our analysis. 

In addition to these galaxies, we have included 31 out of 
the 38 face-on spirals analyzed in ISanchez et al . (2012b). This 
sample comprises face-on spiral galaxies extracted from the 
CALIFA feasibility studies dMarmol-Oueralto et al.ll201 ll) . All 
these galaxies were observed with PPAK, using a setup similar 
to the one used for CALIFA, with slightly different spectral res- 
olutions and wavelength coverage (see Marmol-Queralto et al. 
2011, for details). They cover the same redshift range, optical 
size and the extent covered by the CALIFA, too. On the other 
hand, t he remaining 7 spirals were extracted from the PINGS 
survey (Rosales-Ortega et alJl2010l) . These galaxies are located 
at lower redshift, being larger in projected sizes, and in some 
cases they have not been completely covered by the IFS data. 
Therefore, we prefer to excluded them from this analysis. 



3. Analysis and Results 

3. 1 . Derivation of the analyzed parameters 

The three main quantities to derive for each galaxy are (i) the 
integrated stellar mass, (ii) the characteristic oxygen abundance 
and (iii) the integrated SFR. Thes e parameters were de rived fol- 
lowing the procedure described in Sanchez et al. (2012b), devel- 
oped for a similar set of data. In the following we briefly outline 
the most important aspects of their derivation. 

3.1.1. Stellar Masses 

We derive the stellar masses using the integrated V-band magni- 
tudes and B - V colors listed in Table [T] an d the average mass- 
luminosity ratio (M/L) described in Be ll & de Jon g (2001). The 
B and V-band magnitudes were derived using the Petrosian mag- 
nitudes provided the SDSS ph otometric ca t alog EL and trans- 
formed using the equations by I Jester et all (120051) . The errors 
provided by the SDSS catalog were propagated to derive the 
corresponding errors for the B and V-band magnitudes. The ex- 
pected systematic error of the transformation, of ~0.02 mag, was 
not taken into account since it will affect equally to all the listed 
values. . A direct comparison with our own derived magnitudes, 
based on a growth-curve analysis and a detail subtraction of 
the local background of the SDSS images shows that the errors 
could be ad maximum ~0.07 mag (Walcher et al., in prep.) The 
integrated luminosities and colors were corrected for the effect of 
dust attenuation, where the internal extinction was derived from 
the multi-SSP analysi s of the stellar continuum summarized in 
San chez et alj d2012al) , prior to derive the stellar masses. 

This procedure to derive the stellar masses is robust and 
straight-forward, and it only requires a few assumptions on 
the properties of the galaxies. We adopt the iBell & de Jongl 
(2001 ) M/L ratios, instead of the more recent ones provided 
by Bell et alj d2003l) for the S PSS colors, fo r consist ency with 
the derivation p erformed in ISanchez et al.l (l2012bb . Indeed, 
Bell et al. (2003) show that both M/L ratios are consistent with 
each other. 

In any case, we perform a sanity check to test the robustness 
of the stellar masses derived. For doing so, we derive the masses 
using four additional methods: (i) we use the full spectral energy 
distrubution (SED) provided by the SDSS photometry, and we 
derive the masses using the stellar populatio n analysis provided 
by the PARADISE code dWalcheret all201ll) . This derivation is 
less affected by possible color biases, since it includes the whole 
optical SEDs. However, the dust attenuation affecting the stellar 
populations was not taken into account; (ii) we derive the stellar 
masses using the K-band photometry extracted from the 2MASS 
photometric catalog Q applying a global aperture correction, 
and adopting the M/L ratio described by Longhetti & Saraccol 
(2009). This derivation is clearly the one that is less affected 
by any possible color biases, the effect of young stellar popu- 
lations and dust effects, a priori. However, it presents its own 
problems, in particlar the uncertainties in the M/L ratio related 
to the unknown contribution of the TP-AGB stars in this wave- 

I || ■ 

length regime (e.g. Walcher et al. 2011). Only 128 objects from 
the ones considered here have photometry in the 2MASS cata- 
log; (iii) w e colle cted the stellar masses recently published by 
iPerez et all (120131) . extracted from the CALIFA datacubes, us- 
ing pycAsso (Cid Fernandes et al., submitted), the 3D imple- 
mentation of the Starlight code (ICid Fernandes et al1l2005l) . for 
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those galaxies in common (75 objects). These masses are the 
most self-consistent with the other parameters presented in this 
article, but they rely on a less straight-forward stellar popula- 
tion decomposition. The analysis takes into account the dust at- 
tenuation; and (iv) we derived the stellar masses provided by 
the stellar population decomposition method used to remove the 
underlying stellar population pixel-by-pixel, using FIT3D. This 
method is, in principle, similar to the previous one. However, 
FIT3D uses a much more reduced library of stellar populations, 
and therefore, it provides, a priory, a less accurate M/L ratios and 
stellar masses. 

Figure [2] illustrates the results of this analysis. For each 
panel, it shows the comparison betw een the stellar mass de- 
rived using the Bell & d e Jond (|2001) relation, and those ones 
described before. In all the cases we find a very good agreement, 
close to a one-to-one relation in most cases. The best agreement 
is with the stellar masses derived using only photometric infor- 
mation, either the full SDSS SED distribution or the 2MASS K- 
band data, with a standard deviation of the difference of cr ~0. 15 
dex. The masses derived using only the SDSS photometry are 
slightly lower, in particular in the low mass regime, which is ex- 
pected since they are not corrected by dust attenuation. For the 
stellar masses derived using the CALIFA spectroscopic informa- 
tion, the best agreement is with those provided by pycAsso. The 
standard deviation of the difference is cr -0.14 dex, which is 
consistent with an spectrophotometric accuracy of the CALIFA 
datacubes of -10% with respect to the SDSS one. Finally, the 
largest dispersion is found when comparing the adopted stellar 
masses with those derived from the multi-SSP analysis provided 
by FIT3D (cr ~0.24 dex), as expected. 

Similar results are found when comparing the masses de- 
rived with the total stellar masses listed in the MPA/JHU catalog 
Kauffmann et al. (2003b). Therefore, if there is a systematic ef- 
fect in the derivation of the stellar masses, this is at maximum of 
the order of the typical error estimated based on the propagation 
of the photometric errors (~0. 15 dex). 

3.1.2. Star Formation Rate 

The SFR was derived for each galaxy based on the proper- 
ties of the ionized gas emission. By construction, CALIFA IFS 
data sample most of the optical extent of galaxies, which al- 
lows to derive the total SFR minimizing any aperture effect. 
The SFR was derived from the integrated Hff luminosity. For 
doing so, we first perform a spectroscopic decomposition anal- 
ysis between the underlying stellar population and the ionized 
gas emission lines described using FIT3 D (Sanchez et al. 2007b; 
iRosales-Ortega et al.l l20TTt ISanchez et al] |2012a|bl). The anal- 
ysis w as performed spaxel-by-spaxel, following ISanchez et al.l 
(l201ll) . Each spectrum in the datacube is then decontaminated 
by the underlying stellar continuum using the multi-SSP model 
derived . Then, we selected those spaxels whose ionization is 
dominated by star- f ormation based on emission line ratios (e.g. 
Sanc hez et al.l l2012b), i.e., the Hn regions described latter. The 
areas dominated by diffuse gas emission have bee n excluded on 
purpo se, since their nature is still unclear (e.g. Thilke ret al.l 
2002, , Iglesias-Paramo et al., in prep.). This does not affect the 
results, since diffuse gas comprises less than 5% of the total in- 
tegrated flux of Ha for our galaxies. We repeated all the cal- 
culations using the SFR integrating over all the FoV, instead of 
just only the H 11 regions without any significant difference. The 
observed He intensities were corrected for reddening using the 
Balme r decrement (Ha/HP ) according to the reddening func- 
tion of lCardeUi et all dl989b . assuming R = A V /E(B - V) = 3.1. 



The theore tical value for the intrinsic Balmer line ratios were 
taken from Osterbrock (1989), assuming case B recombination 
(optically thick in all the Lyman lines), an electron density of 
n e - 100 cirT 3 and an electron temperature T e - 10 4 K, 
which are reasonable assumptions for typical star-forming re- 
gions. The values of the SFR were derived for each galaxy adopt- 
ing the classical relations by Kennicutt ( 1998), based on the dust- 
corrected Ha luminosities. 

The corresponding errors for the SFR were derived, for each 
galaxy, by propagating through the equations the errors pro- 
vided by the fitting procedure for the deriv ation of Ha and HS 
emission lines, spaxel-to-spaxel (FIT3D, Sanchez et alj 1201 U 
2012b). These errors do not include the possible zero-point cal- 
ibration errors in the CALIFA datacubes, which are estimated 
to be lower than a ~ 15% through all the wavelength range 
(lHusemannetalJl2013h . 

3.1.3. Oxygen Abundance 

The characte r istic oxygen abundance was defined by 
Zarit skv et al.l (1 19941) as the abundance at 0.4p25, being 
representative of the average value across the galaxy, as we 
indicated before. This distance corresponds basically to one 
effective radius for gala x ies in t he local universe. We demon- 
strated in San chez et al.l (l2012bl) . that the effective radius is a 
convenient parameter to normalize the abundance gradients 
in galaxies at different redshifts. We will address the study of 
abundance gradients and the relation between the integrated 
and characteristic oxygen abundance in a forthcoming article 
(Sanchez et al. in prep.). Thus, we re-define the characteristic 
or effective oxygen abundance as the corresponding value at the 
effective radius. 

To derive this parameter we analyze the radial gradients of 
the oxygen abundance i n all the galaxi e s of our sample, using the 
procedures detailed in Sanc hez et al.1 (1201 2bl) . In a companion 
article we will describe the procedure for this particular d ataset, 
summarizing the main properties of the ionized regions across 
the galaxies (Sanchez, in prep.). We present here just a brief 
summary of the different steps included in the overall process: 
(i) First we create a narrow -band image of 120A width, centered 
at the wavelength of Ha shifted at the redshift of the targets. 
The narrow-band image is properly corrected for the contami- 
nation of the adjacent continuum; (ii) The narrow band image 
is used as an input for HIIexplorepQ an autom atic Hn region 
detection code created for this kin d of analysis (ISanchez et al.l 
|20 1 2bt IRosales-Ortega etal]|2012l) . The code provides with a 
segmentation map that identifies each detected ionized region. 
Then, it extracts the integrated spectra corresponding to each 
segmented region. Figure Q] illustrates the process, showing the 
Ha intensity maps and the corresponding segmentations, for two 
objects. A total of 3435 individual ionized regions are detected 
in a total of 137 galaxies from the sample; (iii) Each extracted 
spectrum is then decontaminated by the underlying stellar con- 
tinuum using the multi-SSP fitting routines included in FIT3D, 
as described before; (iv) Then each emission line within the con- 
sidered wavelength range is fitted with a Gaussian function to 
recover the line intensities and ratios; (v) These line ratios are 
used to discriminate between different ionization conditions and 
to derive the oxygen abundances for each particular ionized re- 
gion. Basically we selected those clumpy ionized regions with a 
strong blue underlying stellar population, that contributes at least 
a 20% of the flux in the V-band, to ensure that they are located 
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Fig. 4. . Top-left panel: Distribution of the oxygen abundances at the effective radii as a function of the integrated stellar masses 
for the CALIFA galaxies (113, blue solid circles). For comparison purposes, we show similar values for the galaxie s observ ed in 
the CALIFA feasibility studies (31, red solid squares) included in the Hn regions catalog described in lSanchez et al.l (1201 2b). The 
size of the symbols is proportional to the integrated SFR. The solid line represents the best fitted curve, as described in the text. 
We also included t he distribution of abunda ces as a function of the stell ar masses for the SPSS data corresponding to the redshift 
ranges sampled bv lLara-Lopez et all Q2010), dashed black contours, and Man nucci et al.l d2010l) . solid grey contours. In both cases, 
the first contour encircles the 95% of the objects, and 20% less objects in each consecutive contour. Top- right panel: Distribution 
of the integrated SFR as a function of the stellar masses. The size of the symbols is proportional to the oxygen abundance shown 
in previous panel. The solid line represents the best fitted linear regression to the CALIFA data. Lik e in the previos panel th e 
contours represent the distri bution of both parame ters corresponding to the redshift ranges sampled bv lLara-Lopez et a l. (2010), 
dashed black contours, and Man nucci et al.l (1201 Oh . solid grey contours Bottom-left panel: Same distribution shown in the top-left 
panel, for the CALIFA and feasibility-studies objects. The colors represent the logarithm of the integrated SFR for each galaxy. 
Bottom-right Panel: Distribution of the differential oxygen abundances at the effective radii, once the dependency with the stellar 
mass has been subtracted, as a function of the integrated SFR for the CALIFA galaxies (113, blue solid circles). The size of the 
symbols is proportional to the oxygen abundances shown in Fig. [4] The histogram shows the same distribution of differential oxygen 
abundances. The solid line represents a Gaussian function with the same central value (0.01 dex) and standard deviation (0.07 dex) 
as the represented distribution, scaled to match the histogram. 



below the lKewlev & Ellison d2008l) demarcation l ine in the clas- 
sical J OIIIl/H/J vs. [NII]/Ha diagnostic diagram (Baldwin et al.l 
1981). A total of 2846 Hn regions are then extracted, for a total 
of 134 galaxies; (vi) Finally, in combination with a morpholog- 
ical analysis of the galaxy, we can recover the abundance gra- 
dients for each particular galaxy, and derive the corresponding 
value at the effective radius. For doing so, we fit the radial distri- 
bution of the oxygen abundance between 0.2 and 2. 1 effective ra- 
dius with a linear regression, and uses the derived zero-point and 
slope to determine the oxygen abundance at one effective radius. 
This procedure increses the accuracy of the derived abundance, 
since it uses the full radial distribution of individual estimations. 
We used only (1) those regions for which the oxygen abundance 



was derived with a nominal error lower than <0.2 dex, which 
corresponds to emission lines with a S/N>5 in all the considered 
emission lines ([OII]3727, HB, [OIII]5007, Ha and [NII]6583), 
and (2) those galaxies with at least 3 Hn regions covering at 
least between 0.3 and 2.1 effective radii. This comprises a total 
of 2061 H ii regions and associations, distributed in 1 13 galaxies. 
This sample comprises galaxies of any type, mostly spirals (both 
early and late type), with and without bars, and with different in- 
clinations. We note that although the detectability of H n regions 
is affected by inclination, i.e., less number is accesible and de- 
tected, we do not notice any further difference in their properties 
(abundance, radial distribution, luminosity...). 
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The oxygen abundance was derived using different strong- 
line indicators, like the R23-fit derived by T04, the N2 and 
03N2 calibrators by Pettini & Pagel (2004), and the recently 
counterpart-method bv lPilvuginet alJ J2012I) . Interestingly, de- 
spite the differences in the derive d estimates, all values show 
a trend with the stellar mass (e.g. iKewlev & Ellisoni 120081) . In 
fact, there is a clear correspondence not only with the calibrators 
(e.g. Lopez-Sanchez et al. 20 12]), but also betwe en the line ra- 
tios used to derive them (e.g. Sanchez et al.ll2012bb . Therefore, 
for simplicity, we will show here only the results based on the 
03N2-calibrator, which is one of the most straight-forward one. 

The results of this analysis are listed in Table Q] including 
the final number of H n regions considered for each galaxy and 
the derived parameters, i.e., the integrated stellar mass and SFR, 
and characteristic oxygen abundance , together with their cor- 
responding errors. In addition to these characteristic values for 
each individual galaxy, we derive similar parameters for each 
individual H n region: the mass surface density (£) and the sur- 
face SFR density (/ySFR). These quantities were derived using 
the same procedures described in Sec. I3.1.1l and 13.1.21 but using 
as observables (i.e., the photometric values and the Ha fluxes), 
those ones corresponding to each particular H n region. Finally, 
the derived values were divided by the corresponding encircled 
area. 



3.2. Mass Metallicity Relation 

We first expl ore the M-Z relation following T04 and 
Mannucc i et alJ d2010l) . Figure |4] shows the distribution of the 
characteristic oxygen abundances as a function of the stellar 
masses for the galaxies of our sample (blue circles). In addi- 
tion, we have i ncluded those values co rresponding to the galax- 
ies analyzed in San chez et alJ (1201 2bl) from the CALIFA feasi- 
bility studies (red squares). As expected, there is a clear trend 
betw een both quantities, foll owing basically the shape described 
by IKewlev & Ellisoni d2008l) for the considered abundance cal- 
ibrator. The correlation coefficient between both quantities is 
r=0.769, indicating that there is a positive trend. From the for- 
mal point of view, this coefficient should be used only for linear 
relations. However, it is still valid as an indication of a trend 
between two parameters. 

For comparison purposes, we include the density distribu- 
tion of galaxies within the M-Z plane for t he SPSS galaxies 
corres ponding to the redshift range studied by Lara-Lopez et al] 
d2010t) (dashed black contours) and Mannu cci et al.l J2010) 
(solid grey contours). The first contour encircles ~95% of the 
total number of galaxies, with a 20% decrease in each con- 
secutive one. We have calculated the oxygen abundances for 
these galaxies using the 03N2 indicator, instead of adopting 
published metallicities which will force us to adopt a correc- 
tion. For doing so, the stellar masses and emission line ratios 
have been directly extracted from the Max-Planck-Institute for 
Astrophysics-John Hopkins University (MPA-JHU) emission- 
line analysis databasqj. The contours show that the CALIFA 
data cover a similar range of masses and metallicities similar 
to both previous studies. However, it is important to note here 
that Lara-Lopez et al. (2010) samples a lowe r mass and higher 
metallicity range than Mann ucci et al] (1201 Ol) . which is a clear 
effect of the redshift range covered by both studies and the SDSS 
selection function. 

Instead of the classical polynomial function adopted by most 
of previous studies for this relation (e.g., T04. IKewlev & Ellison! 



20081: 



Man nucci et al] 

2012) 



120101: iRosales-Ortega et alj 120121: 



Hughes et al. 2012), we prefer to adopt an asymptotic function 
to describe the relation between both parameters following 
Moustaka s~et~ai1 J201ll) : 



y = a + b(x - c) exp[-(x - c)] 



(1) 



Where y = 12+log(0/H), and x is the logarithm of the stellar 
mass in units of 10 8 M G . This functional form tries to describe 
the distribution of parameters defining a maximum abundance 
for large masses (a). Although the adopted formula is slightly 
different for computatio nal convenience than the one adopted by 
Moustak as et al.l d201ll) . the main motivations and interpretation 
are the same. Obviously, due to the sample construction, we do 
not have enough data in the low-mass/low-metallicity range to 
explore this second range. 

The derived parameters for the CALIFA data are 
a=8.74±0.01, Z?=0.0 18 +0.007 and c=3.5+0.3. These two 
latter parameters govern how the metallicity actually increases 
with mass, while the first one is the asymptotic oxygen abun- 
dance for large masses at the effective radius of the galaxies. 
The derived values cannot be extrapolated to masses lower than 
10 9 Mo, not cov ered by the CALIFA m other sample. A recent 
compilation by [Kudritzki et al. (2012) of the characteristic 
metallicities on a set of local universe galaxies indicate that 
while at high mass they reach an asymptotic value, at low 
masses the rela tion is almost linea r, consistent w ith early 
derivations (e.g. ISkillman et al.lll989t ISkillmanl[T992h . With 
the adopted formula it is possible to reproduce this behavior, 
although the parameters should be adjusted. 

From a theoretical point of view the existence of an asymp- 
totic value is predicted, as proceeding from the production of el- 
ements in stars and a given IMF, that is, the stellar true yield for a 
single stellar population which will be reached when all stars die. 
An abundance at the effective radius of ~8.74 dex corresponds to 
a maximum abundance of ~8.84, in the central re gions, in agree- 
ment w ith theoretical expectations. For example. iMolla & Diazl 
(2005) found a saturation level in 12+log(O/H)=8.80. They used 
a grid of chemical evolution models for spiral and irregular 
galaxies, the stellar yields by Wooslev & Weaver ( 1995) for oxy- 
gen and the IMF bv lFerrini et al.l (j.1992). If a different combina- 
tion of IMF and stellar yield set is used, differences of a factor 
two may be found. An average value of 8.80±0. 15 would be the 
maximum oxygen abundance predicted theoretically, which is is 
also in agreement with more recent estimations Pilvugin et al. 
(120071) . 

Since abundance is a relative parameter (amount of oxygen 
relative to the amount of hydrogen), this value can also be inter- 
preted as the maximum abundance of the universe at a certain 
redshift. As shown by Moustaka s et al] (1201 ll) . this parameter 



evolves with redshift at about 



Mog(OIH) 



-0.2. The functional 
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form describes well the considered data, as it can be appreci- 
ated in Fig. |4] The dispersion of the abundance values along 
this curve is <x Alog(0 /H) =0.07 dex, much smalle r than previous 
reported values by (e.g., ~0.1 dex, T04 Mann ucci et al.1 12010: 
Hugh es et al.ll2012l) . It is important to note here that this disper- 
sion is not much larger than the estimated errors for the charac- 
teristic abundances (ei og (o/H) ~0.06 dex). Although the disper- 
sion around the correlation depends on the actual form adopted 
for the M-Z relation, the derived value is always low. If in- 
stead of the considered functional form we a dopt a 3 or 4 order 
polynomial function, as the one considered bv lKewlev & Ellisoni 
(2008), the dispersion is just 0.084 dex, smaller than the value 
reported in that publication for SDSS data using the same cal- 
ibrator. In any case, our adopted functional form has a more 
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Fig. 5. . Top-left panel: Distribution of the oxygen abundances for the 3000 individual H n regions extracted from the CALIFA data, 
as a function of the surface mass density, represented with a density contour-plot. The first contour encircles 95 % of the total number 
of H ii regions, with ~20% less in each consecutive contour. The solid-yellow points represent the average oxygen abundances, with 
their corresponding standard deviations indicated as error bars, for consecutive bins of 0. 1 dex in surface mass density. The dot- 
dashed line represents the best fitted curve for these points. Top-right panel: Distribution of the SFR surface density as a function of 
the stellar mass density of each individual H n region. The contour plots represent the density of points in the distribution of oxygen 
abundances as a function of the surface SFR densities, with the same encircled numbers as in the previous panel. Bottom-left panel: 
Same distribution shown in the top-left panel, for the CALIFA and feasibility-studies objects. The colors represent the logarithm of 
the integrated surface density of the SFR for each H n region. Bottom- right panel: Distribution of the differential oxygen abundance, 
once the dependency with the surface mass density has been subtracted, for individual H n regions. The contour plots represent 
the density of points in the distribution of differential oxygen abundances, once subtracted the dependence on the surface masses 
densities, as a function of the surface SFR densities, with the same encircled numbers as in the previous panel. The solid-yellow 
points represent the average values for the abundance offsets, with their corresponding standard deviations indicated as error bars, 
for consecutive bins of 0. 1 dex in surface SFR density. 



straighforward physical interpretation than the pure polynomial 
functions adopted in the literature. 

The results do not vary significantly when using the oxygen 
abundance derived from the integrated spectra across the covered 
FoV. However, both the dispersion around the derived M-Z rela- 
tion ( ~0. 1 dex ) and the estimated errors for the abundances are 
slightly larger ( ~0.07 dex). This is expected since for the inte- 
grated spectra we coadd regions with different ionization sources 
(e.g. Sanchez et al. 2012b), on one hand, and we have a single 
estimation of the abundance per object, not a value derived from 
the analysis of a well defined gradient (which reduces the error 
in the derivation of the oxygen abundance). Figure |3]illustrates 
this effect, showing the distribution of the oxygen abundances 
estimated for each galaxy at the effective radius, using the radial 
gradients, with the oxygen abundances derived directly from the 
spectra integrated over the entire field-of-view. As expected both 



parameters shows a tight correlation (cr ~0.06 dex). The errors 
estimated for the integrated abundance are larger, and they are 
directly transfereed to the dispersion along the one-to-one rela- 
tion, and thus, to the corresponding M-Z distribution. 

Once this functional form had ben derived, we determined 
the offset and d ispersion for th e data co rresponding to the galax- 
ies analyzed in ISanchez et al] (1201 2b). Based also on IFS, and 
covering the full optical extent of the galaxies, the only pos- 
sible differences would be either the details of the sample se- 
lection and/or the redshift range. The galaxies from the fea- 
sibility studies considered here are all face-on spiral galaxies, 
with an average redshift similar to that of the CALIFA data 
(z ~0.016), and only slightly larger projected sizes (red squares 
in Fig.B). They present a minor offset with respect to the derived 
M-Z relation (Alog(0/H) - -0.003 dex) and a slightly larger 
dispersion (o~Ahg(.0/H) =0.084 dex). Considering that the spec- 
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trophotometric calibration of the data is not that accurate (~20%, 
Marmol-Oueralto et al. 201 1), the slightly larger redshift range 
covered by this sample (Sanchez et al. 2012b), and that the sam- 
ple is smaller (3 1 objects), the difference is not significant. If we 
consider both this dataset and the CALIFA one described before, 
the derived dispersion is just cr^i og ^oiH) =0.07. 

3.3. Dependence of the M-Z relation on the SFR 

The main goal of this study is to explore whether there is a de- 
pendence of the M-Z relation on the SFR. Before doing that, it 
is important to understand the possible additional dependence of 
the SFR on other properties of the galaxies, in particular on the 
stellar mass. The SFR of blue/star-forming galaxies has a well 
known relation with the stellar mass which has been described 
not only in the Local Universe (e.g. Brinchmann et al. 2004), 
but also at different redshifts (e.g. Elb az et aDl2007 , and refer- 



ences there in). Its nature is still under debate (e.g. lDave] [2008: 
IVulcani et al.ll2010l) . but it is most probably related to the ability 
of a certain galaxy to acquire neutral gas, together with the well 
known Schmidt-Kennicutt law (Kennicutt 1998). Regardless of 
its nature, this relation, together with the M-Z relation may in- 
duce a secondary correlation between the SFR and the oxygen 
abundance. 

Figure [4] top-right panel, shows the distribution of the SFR 
as a function of the stellar mass for both the CALIFA galax- 
ies analyzed here and the galaxies extracted from the feasibility 
studies. As expected, both parameters show a clear correlation, 
with a strong correlation coefficient (r =0.77). The slope of this 
relation, for our galaxies, is a =0.63, only slightly smaller than 
the val ue reported for the Local Universe based on SDSS data 
(~0.77. lElbaz et afll2007T) . This value is somewhat s maller than 
the va lue near to one expected from the simulations (lElbaz et alJ 
2007). The limited number of low-SFR objects in our sample, 
may explain the difference in the slope . Nevertheless it is clear 
that we reproduce the expected correlation . In order to inves- 
tigate if the selection criteria of the CALIFA galaxies, or the 
limited number of observed galaxies, have induced any spuri- 
ous correlation that may affect the distribution of both quantities 
and/or the described trend, we have compared these distribution 
with those of the complete SPSS sa mpl e for the redshift range s 
covered bv lLara-Lopez et all (l2010h and lMannucci et all (l2010h . 
For doing so, we uses the MPA-JHU catalog described in the 
previous section, extracting the same stellar masses and the aper- 
ture corrected SFR for both redshift ranges. The distribution of 
both parameters were represented in Fig. |4] top-right panel, us- 
ing similar density contours as the ones shown in the top-left 
panel. It is clear that the ranges of both quantities, and the distri- 
bution along the Mass-SFR plane, agree between the three sam- 
ples. In particular, the agreement is remarkable between our an- 
alyzed CALIFA sample and the one studied bv lLara-Lopez et all 
d2010l) . which is the one at more similar redshift range. Indeed, 
the described correlation agrees completely with the SDSS data 
corresponding to these redshift range. Thus, the selection criteria 
does not seem to affect neither the range of parameters nor the 
correlation between both quantities. 

FigureH] bottom-left panel, shows the M-Z relation, but with 
a different color-coding indicating the integrated SFR of each 
particular galaxy. There is a trend of the SFR with the metallic- 
ity, which follows the trend of the SFR with mass. Thus, more 
massive galaxies have both higher SFR and are more metal rich. 
In order to determine which is the dominant correlation, we de- 
rive the correlation coefficient between the SFR and the oxygen 
abundance, obtaining a value of r = 0.54. This indicates that, 



while the strength of the M-Z and the SFR-Mass relations are 
indistinguishable statistically, the SFR-abundance one is weaker, 
and seems to be induced by the other two. 

We do not find any trend between the SFR and the oxy- 
gen abundance for a give n mass, as the one d escribed by 
dLara-Lopez et al.l2010l) andlMannucci et al] (l2010l) . For a given 
mass, the galaxies with stronger star-formation do not seem to 
have a lower metal content. In order to investigate it further, 
we derived the differential oxygen abundance with respect to the 
considered M-Z relation, to search for any correlation with the 
SFR. 

Figure [4] bottom-right panel, shows the distribution of the 
residual of the oxygen abundance, once the M-Z relation has 
been subtracted, as a function of the SFR. The data is consistent 
with being randomly distributed around the zero value, with no 
correlation between both parameters (r = 0.08). The histogram 
of the differential oxygen abundances shows a distribution com- 
patible with a Gaussian function. We recall here that the disper- 
sion around the M-Z relation was just ~0.07 dex, a markedly 
small error taking into account the nominal errors of the param- 
eters involved. In summary, we cannot confirm the existence of 
a dependency of the M-Z relation on the SFR, with the current 
CALIFA data. 

A possible concern regarding this result may be that the lim- 
ited number of objects in the current CALIFA sample may ham- 
per the detection of a significant correlation. We perform a sim- 
ple Monte-Carlo simulation to test this possibility. For doing 
so, we assume that the M-Z -SFR distribution follows the an- 
alytical formulae described by Mannucci et al. (2010). Then, we 
extract a random subsample of 100 objects covering a similar 
range of the parameters covered by the CALIFA sample, and 
we add a gaussian n oise of 0.05 dex (i.e., the value reported by 
iMannucci et al.ll201 0). For this subsample, we subtracted the M- 
Z relation predicted by the same authors, and analyze the possi- 
ble correlation between the residual of the abundance and the 
star- formation rate. After 1000 realizations we found that in all 
the cases we found a clear correlation (r ~ 0.7), with an slope 
compatible with the one predicted by the FMR relation. Thus, 
the number of objects does not affect the result. 

3.4. Dependence of the L-Z relation on the SFR 



In lRosales-Ortega et alJ (J2012) we presented a new scaling rela- 
tion between the oxygen abundance and the surface mass density 
for individual H n regions (E-Z relation)Q.This relation is a local 
version of the global M-Z relation, and indeed, it is possible to 
recover this latte r one from the former, just b y integrating over a 
certain aperture (Rosales-Ort ega et al.l l2012). We proposed that 
this relation is the real fundamental one and the M-Z relation 
being a natural consequence. Following this reasoning, it is im- 
portant to test if there is a secondary relation between the local- 
M-Z relation and the SFR. An additional advantage of this ex- 
ploration is the increase of statistical significance of the results, 
due to the much larger number of sampled objects here (~3000 
Hn regions). 

Figure [5] top-left panel, shows the E-Z relation for the 
H n regions of the CALIFA galaxies explor ed so far. As a first re- 
sult, w e reproduce the relation described in lRosales-Ortega et alJ 
(2012), using a completely different dataset (the catalog of 
Hn regions extracted from 38 face-on spirals, using IFS data, 



An 



online 



annimated 



of this 



rela- 



tion for the current analysed data is available at 
http : //www . youtube . com/watch?v=F7 JCX- d7uPY 
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Fig. 6. Left panel: Same distribution of the characteristic oxygen abundance as a function of the stellar mass as shown in the 
bottom panel of Fig. [4] Right panel: Same distribution of the oxygen abundances as a function of the stellar surface density for the 
individual H n regions, as shown in the bottom panel of Fig. In both cases, the color code indicates the logarithm of the specific 
star formation rate (sSFR). 



by ISanchez et al.|[2012bl) . Both quantities have a strong corre- 
lation, with a correlation coefficient of r=0.98. As in the case 
of the global M-Z relation, we adopted the asymptotic equation 
Q] In this case, the x parameter represents the logarithm of the 
stellar mass surface density in units of M pc~ 2 . The a, b and 
c parameters have the same interpretation as those derived for 
the global relation, and we derive similar values when fitting the 
data: a =8.86+0.01, fc=0.97±0.07 and c=2.80+0.21, reinforcing 
our hypothesis that both relations are strongly linked. In general, 
the global-A4-Z relation has the same shape as the local one but 
shifted arbitrarily in mass (to match them to the corresponding 
mass densities). 

Like in the case of the integrated properties of the galax- 
ies, there is a tight correlation between the stellar surface mass 
density and the surface star-formation density (Fig. top-right 
panel). The correlation is a clear correlation (r = 0.66), and it 
is well described by a simple linear regression. The surface star- 
formation density scales with the surface mass density with a 
slope of a=0.66+0.18, once more, below the expected value de- 
rived in semi-analytical simulations for the Mass-SFR relation. 
Note that a comparison with a simulation taking into account 
the described local relationship is not feasible at present, since 
no such simulation exists to our knowledge. 

This relation indicates that the HII regions with stronger star- 
formation are located in more massive areas of the galaxies (i.e., 
towards the center). We have explored if our result is affected 
by the presence of young stars that biases our mass derivation. 
However, we found that our surface-brightness and colors (the 
two parameters used to derive the masses), are not strongly cor- 
related with either the age of the underlying stellar population or 
the fraction of young stars. Evenmore, the SFR-density presents 
a positive weak trend with the B-V color, not the other way 
around. 

Like in the global case, these two relations between the oxy- 
gen abundance and the surface star- formation rate with the sur- 
face mass density, may induce a third relation between the two 
former parameters. In fact, there is a trend between both param- 
eters, much weaker than the other two described here (r - 0.59), 
as expected from an induced correlation. Figure |5] bottom left 
panel illustrates this dependence, showing the local- M-Z rela- 
tion, color-coded by the mean surface star- formation rate per bin 
of oxygen abundance and surface mass density. There is clear 



gradient in the surface SFR which mostly follows the local- 
M-Z relation. In this figure we can see a de viation from this 
global tren d in agreement with th e results by lLara-Lopez et al.l 
(2010) and Mannucci et al. (2010), in the sense that there are a 
few H ii regions with low metal content and high star-formation 
rate for their corresponding mass surface density (bottom right 
edge of the figure). However, all together they represent less than 
a 2-3% of the total sampled regions, and therefore cannot ex- 
plain a global effect as the one described in those publications. 
Furthermore, the stochastic fluctuations around the average val- 
ues at this low number statistics may clearly produce this effect, 
since there are other H n regions in the same location of the 
diagram with much lower surface SFR. 



Following the same procedure outlined in Sec. 13.31 we re- 
moved the dependence of the oxygen abundance on the surface 
mass density by subtracting the derived relation between both 
parameters. Then, we explore possible dependencies on the sur- 
face density SFR. Figure [5] bottom-right panel, illustrates the 
result of this analysis. It shows the distribution of the differential 
oxygen abundance along the surface star-formation rate for the 
individual Hn regions. No correlation is found between both 
for the individual values (r = 0.17). However, the dispersion 
around the zero-value is larger than that reported for the similar 
residual distribution for the global relation (cn oca i =0.11 dex). 
The distribution seems to show some tail towards higher abun- 
dances for weaker SFRs and lower for stronger ones. Thus, in 
order to explore possible dependencies even more, we derive 
average differential oxygen abundance for equal ranges of sur- 
face SFR of 0.2 dex width, and determined the correlation coef- 
ficient. In this case we found a clearer trend (r - 0.67), in the 
sense that H n regions with higher surface SFR shows slightly 
lower abundances for similar masses. However, the slope of the 
derived regression a = -0.0155±0.217 is compatible with no 

trend at all, and it is far below the reported value of a 0.3 for 

the dep endence of the M-Z relation and the SFR Mannucci et al.l 

(1201 Oh or a 0.4 bv lLara-Lopez et all d2010t) . Applying any 

of these correlations (i.e., adding fflog(SFR) to the mass) does 
not reduce the dispersion found around the zero value at all 
(o~i cai,SFR =0.108 dex). In summary, we do not find a correla- 
tion between the local-A4-Z relation with the SFR, in the way 
described by other authors for the global one. Our results for the 
local and global M-Z relations are therefore fully consistent. 
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Fig. 7. Distribution of the specific star formation as a function of 
the absolute equivalent width of Ha, for the individual H n re- 
gions. The solid yellow circles indicate the location of the aver- 
age values in bins of 0.2 dex of the equivalent width of Ha. The 
solid line shows the best fitted linear regression line 



3.5. Dependence of the M-Z relation on the sSFR 



Following Mannucc i et al.l (l2010h . we have explored the possible 
dependence of the M-Z relation on the specific star-formation 
rate (sSFR), a measurement of the current star-formation activ- 
ity of the galaxy in comparison with the past one. It is known to 
exhibit a smooth dependence on the stellar mass (sSFRocMr 1 ), 
which has been described at different redshifts (e.g. Elbaz et al. 
2007). Ma ny properties of galaxies depend on this parameter. In 
particular, Mannuc ci et al.l (1201 Ot) found that there was an appar- 
ent threshold in the sSFR at about 10~ 10 yr _1 which discrimi- 
nates the effect of the SFR on the M-Z relation, in the sense that 
the depence is stronger for higher values of the sSFR. In this sec- 
tion we explore the possible dependence of the M-Z on sSFR. 

Figure [6] shows both the global and local M-Z relations for 
the galaxies and H n regions discussed in this article. In this case 
the color code indicates the sSFR, determined by dividing the 
SFR and stellar masses described before. In both cases the sSFR 
seems to follow the dependence of the abundance with the mass, 
without a significant deviation from the trend along this relation. 
Thus, we do not identify neither a clear secondary dependence 
on the sSFR, nor a clear threshold at which the global and local 
M-Z relation show a change in the trend or an offset with respect 
to the average relation. 

In the case of the E-Z relation, the data seem to fol- 
low a smooth surface in 3D space, with the three parame- 
ters clearly correlate d . This surface was already discussed by 
iRosales-Ortega et aD (120121) . who, instead of using the sSFR 
used the equivalent width of Ha, a w ell-known proxy for the for- 
mer parameter (e.g. lKennicutJ[l998h . Figure|7]shows the empir- 
ical relation for both parameters based on the individual H n re- 
gions. The correlation is very strong (r - 0.84), and well de- 
scribed with a simple linear regression, with a dispersion 0.26 
dex: 

log(sSFR/yr) = -11.38 ±( ,oi + 1.23 ±0 .o2 • log|EW(Ha)/A| (2) 

Adopting this relation, the described curve is basically the same 
as discussed by Rosales-Ortega et al. (2012), considering the 
reparametrization of one of the axes. We note here that it is re- 
quired to take the absolute value of the equivalent width of Ha, 



which has always negative values for strong emission line re- 
gions. 

4. Discussion and conclusions 

Along this study we have explored the global and local M-Z rela- 
tions and their possible dependence on the SFR and sSFR based 
on the ionized gas properties of the first 150 galaxies observed 
by the CALIFA survey. So far, we have not found any depen- 
dence on the SFR and sSFR different than that induced by the 
well known relation of both parameters with the stellar mass. 
In addition, we explored similar dependencies for the sample 
of ~3000 individual H n regions extracte d from the IFS data. 
We co nfirm the E-Z relation described in iRosales-Ortega et all 
(120121) . However, no secondary dependence is found between 
this relation and the SFR (or sSFR), different than that induced 

by their dependences on the stellar mas_s. 

Our results co ntradict those of Mann ucci et al.l (1201 Ol) and 
lLara-Lopez et al.l (12010). who described a trend for which galax- 
ies with stronger star-formation show lower oxygen abundances 
for the s ame mass range. Si milar results were already pre- 
sented by Hughes et al.l d2012l) . using a sample of galaxies ob- 
served using drift-scan technique s. The main difference be- 
tween both datasets is that while Mannu cci et alj d2010l) and 
lLara-Lopez et al.l d2010l) use the spectroscopic data provided by 
the SPSS (i.e ., singl e aperture, with a strong aperture bias), both 
Hug hes et al.l d2012l) and the results presented here are based 
on integrated (and spatially resolved) properties of the galax- 
ies. Shen et al. (2003) show that the physical size of galaxies 
in the SDSS survey comprises a range of values between a few 
kpcs to a few tens of kpcs. This means that the SDSS fiber cov- 
ers a range between ~0.3 and ~10 tim es the effective radii of 
galaxies at the redshift range covered by iMannucci et al.l (120 10) 
and Lara-Lopez et al. (2010), depending both on their intrinsic 
properties and the involved cosmological distances. This aper 
ture bias affects mostly the derived SFR (depending on the aper- 
ture corretion applied, which are under debate, e.g.lGerssen et al 



1201 2t) and the oxygen abundances (e.g. Sanchez et alJl2012af 
but in some cases it may also affect the mass derivation (if the 
spectroscopic information is used to derive the Mass-to-Light 
ratio). 

In addition, the redshift ranges co vered by these studie s are 
also different. The data presented by Hug hes et al.l (1201 2l) and 
ours are focused on relatively nearby galaxies (D < 120Mpc), 
with a narrow ran ge of distances in both cas es. In contrast, 
SPSS data used by Mann ucci et alJ d2010h and lLara-Lopez et all 
(120101) also differ (0.07 < z < 0.30 and 0.04 < z < 0.1, respec- 
tively), however they both correspond to a much wider redshift 
range than that covered by us. In both cases the aperture biases of 
the SPSS spectra have two effects: (i) They grab a larger portion 
of the galaxies at higher-redshifts, for the same masses, which 
affects their estimations of the abundance (decreasing the value 
at higher redshifts) and the SFR (increasing the value at higher 
redshifts too); (ii) Pue to the reported redshift evolution of the 
M-Z relation, they derive a lower abundance for higher redshift 
galaxies than for lower-redshift ones, for the same mass, even 
in the absence of any aperture bias. This second effect may 
be less importan t, since it would be a t maximum of the order 
of -0.06 dex forlMannucci et al.l (l2010t) . but just -0.02 dex for 
lLara-Lopez et al.l d2010l). The first effect is also stronger for the 
case of the study by IMannucci et all (12010). We can reproduce 
qualitatively the secondary trend with the SFR of the M-Z rela- 
tion by taking our IFS data and simulate different SPSS aper- 
ture spectra corresponding to different projected distances (see 
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Appendix [AJ. However, the derived slope is not reliable due to 
the uncertainties involved (different sample selections, different 
methods to derive the metallicity, unclear evolution of Al-Z re- 
lation with the redshift). Although we cannot unambiguously 
prove that aperture effects are inducing this relation, qualitatively 
there is certainly an agreement. 

As described by iMannucci et al.l (1201 Ol) . oxygen abundance 
is a rather simple parameter, which is governed by three pro- 
cesses: (i) star-formation, which transforms gas into stellar mass, 
increasing the metal content; (ii) gas accretion, or inflows, which 
dilutes the oxygen abundance (assuming accretion of more or 
less pristine gas); and (iii) gas outflows, induced by the star- 
formation. However, none of these processes is simple, and their 
interdependence is unclear. 

Based on the Schmidt-Kennicutt law, the surfac e SFR de- 
pends on the surface mass density, yuSFR oc ~L gas n dKennicutll 
1998), with a slope larger than one (n ~ 1.4 - 1.5). This means 
that the integrated SFR of a galaxy should depend on the gas 
inflow rate (e.g. Mannucci et al. 2010). However, this does not 
necessarily mean a global dilution and a decrease in the aver- 
age oxygen abundance, since inflow is a radial process while 
star-formation is a local one. Gas that falls into the inner regions 
may be polluted by star formation further out if the latter pro- 
cess is fast enough. Therefore, if the typical timescale to recycle 
gas is shorter than the inflow timescale (e.g. Silk 1993), no SFR 
dependence is expected. On the other hand, if the timescales 
for the feedback and global gas recycling are much longer than 
the typical length of a star-forma tion episode (~1 Gyr for L* 
galaxies Ouillen & Bland-Hawthorn 2008), there will be a lose 
of connection between the starformation and the global metal 
enrichment, too. 

On the other hand, outflows induced by strong star- formation 
and the subsequent winds due to SN explosions are known to 
eject metals into the intergalactic med ium, and therefore regu- 
late the metal content of gal axies ("e.g. lLarsonl[l97 4 ; Prieto et all 
120081; iPalacios et"ai] 120051) . It is therefore expected that the 
strength of the outflows, and the fraction of metals lost, depends 
on the SFR. However, the detailed relation is difficult to address, 
since the ability of an outflow to eject metals could be counter- 
balanced by the gravitational field of the galaxy, which makes 
metals to fall back into the disk in a back rain process. The high 
metallicity of the diffuse IGM in galaxy groups and clusters is in- 
direct evidence that at least some fract ion of the metals i n these 
winds is lost to the environment (e.g. Renz ini et alJll993l) . 

Although it is hard to study just using the integrated prop- 
erties of galaxies, the spatial distribution of metals and, in par- 
ticular, the dispersion of the oxygen abundance at a certain dis- 
tance from the center could be used to constrain the strength 
of these (and other) diffusion processes. Our preliminary results, 
Sanc hez et al.l d2012bl) . indicate that this dispersion is rather low. 
Therefore, either metals are really expelled and lost and/or the 
net effect is not very strong. The lack of a clear secondary de- 
pendence of the oxygen abundance, both globally or locally on 
the SFR and the //SFR also suggests that this is not the dominant 
process. 

Fo llowing the formalism presented by IMannucci et al.l 
d2010l) . the fraction of metals ejected due to outflows has to be 
compensated by the metals produced by the star-formation pro- 
cess. If, instead of the reported dependence of the abundance on 
the SFR of ~ -0.32-log(SFR), we consider no dependence at 
all, and using their formulae, the amount of metals lost by this 
process should be Alog(0/H) oc SFR 33 , i.e., much smaller than 
that required to reproduce their results. 



In summ ary our results are c onsistent with those already 
presented in ISanchez et al.1 (1201 2bl) and Rosal eiTOrtega et all 
(2012). In general, the properties of the ionized gas in late-type 
galaxies are consistent with a quiescent evolution, where gas re - 
cycling is faster than other times scales involved (ISilk| [l993). 
This would imply that the galaxies seem to behave locally in a 
similar manner than globally, dominated by a radial mass distri- 
bution following the potential well of the matter, with an inside- 
out growth that is regulated by gas inflow and local down-sizing 
star-formation. Therefore, the dominant parameter that defines 
the amount of metals is the stellar mass, since both parameters 
are the consequence of an (almost) closed-box star-formation 
process. However, we should remark that our results are only 
valid for galaxies with stellar masses higher than ~10 95 M , 
where the CALIFA sample becomes complete. 
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Table 1. Properties of the sample of galaxies. 



Galaxy 


redshift 


V-band (mag) 


B-V 


Type 


Nhii 


log(Mass/MO) 


12+log(0/H) 


log(SFR/MOyr') 


IC5376 


0.01663 


13.84 ± 0.02 


0.89 ± 0.03 


SbA 


2 


10.55 ±0.19 






NGC7819 


0.01730 


13.80 ± 0.02 


0.71 ±0.03 


ScA 


31 


10.41 ±0.17 


8.63 ± 0.07 


0.48 ± 0.01 


UGC00036 


0.02089 


13.72 ±0.01 


0.95 ± 0.01 


SabAB 


5 


10.88 ± 0.20 


8.72 ± 0.04 


0.40 ± 0.01 


NGC0001 


0.01502 


13.11 ±0.01 


0.85 ± 0.01 


SbcA 


31 


10.75 ±0.19 


8.75 ± 0.06 


0.93 ± 0.01 


NGC0036 


0.01993 


13.06 ± 0.01 


0.93 ± 0.01 


SbB 


25 


11.08 ±0.20 


8.71 ± 0.04 


0.76 ± 0.01 


UGC00312 


0.01532 


13.73 ± 0.02 


0.48 ± 0.02 


SdB 





10.06 ±0.14 






NGC0444 


0.01596 


14.34 ± 0.03 


0.63 ± 0.04 


ScdA 


11 


9.98 ±0.16 


8.52 ± 0.06 


0.10 ±0.01 


NGC0477 


0.01989 


14.05 ± 0.02 


0.81 ±0.02 


SbcAB 


24 


10.53 ±0.18 


8.64 ±0.11 


0.75 ± 0.01 


IC 1683 


0.01584 


13.67 ±0.01 


0.88 ± 0.02 


SbAB 


4 


10.55 ±0.19 


8.74 ± 0.03 


0.77 ± 0.01 


NGC0496 


0.02003 


13.71 ±0.02 


0.68 ± 0.02 


ScdA 


30 


10.57 ±0.17 


8.62 ± 0.07 


0.81 ±0.01 


UGC01057 


0.02107 


14.23 ± 0.02 


0.62 ± 0.03 


ScAB 


21 


10.33 ±0.16 


8.55 ±0.15 


0.60 ± 0.01 


NGC0776 


0.01616 


13.00 ± 0.01 


0.90 ± 0.01 


SbB 


34 


10.91 ± 0.20 


8.77 ± 0.07 


0.88 ± 0.01 


UGC01938 


0.02118 


14.23 ± 0.03 


0.72 ± 0.04 


SbcAB 


15 


10.44 ±0.17 


8.60 ± 0.08 


0.57 ± 0.01 


NGC1056 


0.01512 


12.56 ± 0.01 


0.92 ± 0.01 


SaA 


10 


10.06 ± 0.20 


8.62 ± 0.03 


0.18 ±0.01 


NGC1167 


0.01640 


12.26 ± 0.01 


0.98 ± 0.01 


S0A 


1 


11.31 ±0.21 






NGC1349 


0.02191 


13.02 ± 0.02 


0.88 ± 0.02 


E6A 


7 


11.13 ±0.19 


8.71 ± 0.08 


0.35 ± 0.02 


UGC03107 


0.02773 


14.41 ± 0.05 


0.82 ± 0.06 


SbA 





10.71 ±0.19 






UGC03253 


0.01363 


13.30 ±0.01 


0.71 ±0.02 


SbB 


26 


10.38 ±0.17 


8.70 ± 0.06 


0.32 ± 0.01 


NGC2253 


0.01180 


12.65 ± 0.02 


0.79 ± 0.02 


SbcB 


20 


10.59 ±0.18 


8.76 ± 0.03 


0.29 ± 0.01 


NGC2410 


0.01570 


12.94 ± 0.01 


0.86 ± 0.01 


SbAB 


19 


10.83 ±0.19 


8.65 ± 0.05 


0.81 ±0.01 


UGC03944 


0.01289 


13.92 ± 0.03 


0.57 ± 0.03 


SbcAB 


34 


9.88 ±0.16 


8.54 ±0.12 


0.03 ± 0.01 


UGC03995 


0.01601 


12.82 ± 0.01 


0.77 ± 0.01 


SbB 


19 


10.84 ±0.18 


8.67 ± 0.05 


0.51 ±0.01 


NGC2449 


0.01629 


13.12 ±0.01 


0.87 ± 0.01 


SabAB 


9 


10.83 ±0.19 


8.77 ± 0.03 


0.33 ± 0.01 


UGC04132 


0.01718 


13.28 ±0.01 


0.80 ± 0.01 


SbcAB 


11 


10.74 ±0.18 


8.67 ± 0.06 


1.18 ±0.01 


UGC04461 


0.01657 


13.87 ± 0.02 


0.55 ± 0.03 


SbcA 


19 


10.14 ±0.15 


8.48 ±0.10 


0.53 ± 0.01 


IC2487 


0.01427 


13.52 ±0.01 


0.84 ± 0.02 


ScAB 


23 


10.51 ±0.19 


8.65 ± 0.06 


0.44 ± 0.01 


NGC2906 


0.01706 


12.51 ±0.01 


0.88 ± 0.01 


SbcA 


31 


10.34 ±0.19 


8.79 ± 0.05 


0.17 ±0.01 


NGC2916 


0.01229 


12.87 ± 0.01 


0.65 ± 0.01 


SbcA 


50 


10.41 ±0.16 


8.72 ± 0.07 


0.47 ± 0.01 


UGC05108 


0.02692 


14.02 ± 0.02 


0.84 ± 0.02 


SbB 





10.87 ±0.19 






UGC05358 


0.01045 


14.99 ± 0.03 


0.64 ± 0.04 


SdB 


6 


9.36 ±0.16 


8.41 ± 0.05 


-0.43 ± 0.02 


UGC05359 


0.02812 


14.38 ± 0.03 


0.73 ± 0.03 


SbB 





10.65 ±0.17 






UGC05396 


0.01778 


14.14 ±0.02 


0.67 ± 0.03 


SbcAB 


16 


10.23 ±0.17 


8.62 ± 0.06 


0.39 ± 0.01 


NGC3106 


0.02074 


13.03 ±0.01 


0.88 ± 0.01 


SabA 


9 


11.08 ±0.19 


8.71 ± 0.07 


0.43 ± 0.02 


NGC3057 


0.01506 


13.86 ± 0.03 


0.43 ± 0.03 


SdmB 


22 


8.95 ±0.14 


8.34 ± 0.08 


-0.49 ± 0.01 


UGC05498NED01 


0.02090 


14.21 ± 0.02 


0.94 ± 0.03 


SaA 





10.66 ± 0.20 






UGC05598 


0.01875 


14.46 ± 0.03 


0.71 ±0.03 


SbA 


10 


10.19 ±0.17 


8.60 ± 0.05 


0.44 ± 0.01 


NGC3303 


0.02124 


13.98 ±0.01 


0.97 ± 0.02 


SOaAB 





10.84 ± 0.20 






UGC05771 


0.02473 


13.67 ± 0.02 


0.96 ± 0.02 


E6A 





11.08 ±0.20 






NGC3381 


0.01638 


13.03 ±0.01 


0.54 ± 0.01 


SdB 


49 


9.58 ±0.15 


8.64 ± 0.06 


-0.05 ± 0.01 


UGC06036 


0.02168 


13.64 ± 0.01 


0.98 ± 0.01 


SaA 





10.99 ± 0.21 






UGC06312 


0.02103 


13.81 ±0.02 


0.93 ± 0.02 


SabA 


4 


10.87 ± 0.20 


8.67 ± 0.04 


0.49 ± 0.01 


NGC3614 


0.01762 


13.15 ±0.02 


0.77 ± 0.02 


SbcAB 


52 


9.95 ±0.18 


8.71 ± 0.09 


0.04 ± 0.01 


NGC3687 


0.01829 


12.80 ± 0.01 


0.65 ± 0.01 


SbB 


53 


10.08 ±0.17 


8.70 ± 0.05 


-0.09 ± 0.01 


NGC3991 


0.01171 


13.72 ±0.01 


0.46 ± 0.01 


SmA 


1 


9.75 ±0.14 






NGC4003 


0.02179 


13.60 ± 0.01 


0.93 ± 0.02 


SOaB 





10.95 ± 0.20 






UGC07012 


0.01030 


14.23 ± 0.02 


0.47 ± 0.02 


ScdAB 


13 


9.47 ±0.14 


8.49 ± 0.08 


-0.06 ± 0.01 


NGC4047 


0.01125 


12.44 ± 0.01 


0.68 ± 0.01 


SbcA 


47 


10.55 ±0.17 


8.71 ± 0.07 


0.71 ±0.01 


UGC07145 


0.02195 


14.32 ± 0.03 


0.71 ±0.03 


SbcA 


20 


10.40 ±0.17 


8.62 ±0.13 


0.76 ± 0.01 


NGC4185 


0.01284 


12.87 ± 0.01 


0.79 ± 0.02 


SbcAB 


30 


10.58 ±0.18 


8.75 ± 0.06 


0.33 ± 0.01 


NGC4210 


0.01893 


12.92 ± 0.01 


0.72 ± 0.01 


SbB 


56 


10.11 ±0.17 


8.73 ± 0.06 


0.06 ± 0.01 



Basel, University of Cambridge, Case Western Reserve University, University 
of Chicago, Drexel University, Fermilab, the Institute for Advanced Study, 
the Japan Participation Group, Johns Hopkins University, the Joint Institute 
for Nuclear Astrophysics, the Kavli Institute for Particle Astrophysics and 
Cosmology, the Korean Scientist Group, the Chinese Academy of Sciences 
(LAMOST), Los Alamos National Laboratory, the Max-Planck-Institute for 
Astronomy (MPIA), the Max-Planck-Institute for Astrophysics (MPA), New 
Mexico State University, Ohio State University, University of Pittsburgh, 
University of Portsmouth, Princeton University, the United States Naval 
Observatory, and the University of Washington. 

This publication makes use of data products from the Two Micron All 
Sky Survey, which is a joint project of the University of Massachusetts and 
the Infrared Processing and Analysis Center/California Institute of Technology, 



funded by the National Aeronautics and Space Administration and the National 
Science Foundation. 



References 

Baldwin, J. A., Phillips, M. M., & Terlevich, R. 1981, PASP, 93, 5 
Bell, E. F. & de Jong, R. S. 2001, ApJ, 550, 212 

Bell, E. E, Mcintosh, D. H., Katz, N., & Weinberg, M. D. 2003, ApJS, 149, 289 
Boissier, S. & Prantzos, N. 1999, MNRAS, 307, 857 

Brinchmann, J., Chariot, S., White, S. D. M., et al. 2004, MNRAS, 351, 1151 
Brooks, A. M., Governato, E, Booth, C. M., et al. 2007, ApJ, 655, L17 
Calura, E, Pipino, A., Chiappini, C, Matteucci, E, & Maiolino, R. 2009, A&A, 
504, 373 



14 



S. F. Sanchez et al.: The Mass-Metallicity relation explored with CALIFA: 
Table 1. Properties of the sample of galaxies. Continue 



Galaxy 


redshift 


V-band (mag) 


B-V 


Type 


N H ii 


log(Mass/MO) 


12+log(0/H) 


log(SFR/MOyr') 


IC0776 


0.01922 


14.55 ± 0.07 


0.63 ± 0.09 


SdmA 


15 


9.43 ±0.17 


8.33 ± 0.06 


-0.69 ± 0.02 


NGC4470 


0.01776 


12.66 ± 0.01 


0.52 ± 0.01 


ScA 


18 


9.87 ±0.15 


8.59 ± 0.03 


0.20 ± 0.01 


NGC4676A 


0.02242 


14.31 ±0.02 


0.94 ± 0.03 







10.71 ±0.20 






NGC4711 


0.01353 


13.49 ± 0.01 


0.69 ± 0.02 


SbcA 





10.27 ±0.17 






UGC08107 


0.02775 


13.85 ± 0.02 


0.89 ± 0.02 


SaA 





11.02 ±0.19 






NGC4961 


0.01845 


13.39 ± 0.01 


0.48 ± 0.01 


ScdB 


37 


9.63 ±0.14 


8.53 ±0.10 


-0.07 ± 0.02 


NGC5000 


0.01844 


13.55 ± 0.02 


0.74 ± 0.02 


SbcB 


27 


10.61 ±0.18 


8.74 ± 0.05 


0.60 ± 0.01 


UGC08250 


0.01729 


14.93 ± 0.04 


0.63 ± 0.05 


ScA 


6 


9.86 ±0.16 


8.50 ± 0.06 


0.10 ±0.01 


UGC08267 


0.02418 


14.44 ± 0.03 


0.91 ± 0.03 


SbAB 


2 


10.69 ± 0.20 






NGC5016 


0.01848 


12.71 ±0.01 


0.63 ± 0.01 


SbcA 


46 


10.09 ±0.16 


8.75 ± 0.07 


0.18 ±0.01 


NGC5218 


0.01974 


12.74 ± 0.01 


0.90 ± 0.01 


SabB 


12 


10.50 ± 0.20 


8.70 ± 0.05 


0.36 ± 0.01 


UGC08733 


0.01772 


14.29 ± 0.03 


0.68 ± 0.03 


SdmB 


12 


9.37 ±0.17 


8.40 ± 0.09 


-0.47 ± 0.01 


IC0944 


0.02310 


13.45 ± 0.01 


0.92 ± 0.02 


SabA 





11.08 ±0.20 






UGC08778 


0.01065 


13.84 ± 0.01 


0.78 ± 0.02 


SbA 


3 


9.99 ±0.18 


8.67 ± 0.01 


-0.44 ± 0.01 


UGC08781 


0.02518 


13.49 ± 0.01 


0.82 ± 0.02 


SbB 


18 


11.02 ±0.19 


8.68 ± 0.04 


0.61 ± 0.01 


NGC5378 


0.01977 


13.05 ± 0.01 


0.87 ± 0.01 


SbB 


6 


10.33 ±0.19 


8.75 ± 0.03 


-0.43 ± 0.02 


NGC5406 


0.01791 


12.73 ± 0.01 


0.84 ± 0.01 


SbB 


49 


11.02 ±0.19 


8.78 ± 0.06 


0.78 ± 0.01 


UGC09067 


0.03099 


13.96 ± 0.02 


0.66 ± 0.02 


SbcAB 





10.81 ±0.17 






NGC5614 


0.01282 


12.16 ±0.01 


0.95 ± 0.01 


SaA 


8 


11.06 ±0.20 


8.72 ± 0.03 


0.23 ± 0.01 


NGC5633 


0.01766 


12.61 ±0.01 


0.65 ± 0.01 


SbcA 


29 


10.04 ±0.17 


8.77 ± 0.06 


0.43 ± 0.01 


NGC5630 


0.01858 


13.39 ±0.01 


0.46 ± 0.01 


SdmB 


14 


9.60 ±0.14 


8.39 ± 0.05 


0.20 ± 0.01 


UGC09291 


0.01957 


13.58 ± 0.02 


0.61 ± 0.02 


ScdA 


27 


9.81 ±0.16 


8.60 ±0.10 


-0.04 ± 0.01 


NGC5656 


0.01051 


12.48 ± 0.01 


0.70 ± 0.01 


SbA 


31 


10.46 ±0.17 


8.72 ± 0.06 


0.49 ± 0.01 


NGC5682 


0.01724 


14.23 ± 0.02 


0.53 ± 0.03 


ScdB 


5 


9.22 ±0.15 


8.41 ± 0.05 


-0.36 ± 0.01 


NGC5720 


0.02588 


13.68 ±0.01 


0.76 ± 0.02 


SbcB 


6 


10.87 ±0.18 


8.61 ± 0.06 


0.82 ± 0.01 


NGC5732 


0.01241 


13.83 ±0.01 


0.57 ± 0.02 


SbcA 


24 


9.92 ±0.15 


8.63 ±0.12 


0.13 ±0.01 


UGC09476 


0.01065 


13.33 ± 0.01 


0.66 ± 0.02 


SbcA 


46 


10.06 ±0.17 


8.68 ± 0.05 


0.29 ± 0.01 


NGC5735 


0.01244 


13.37 ± 0.02 


0.75 ± 0.02 


SbcB 


52 


10.33 ±0.18 


8.65 ±0.10 


0.30 ± 0.01 


NGC5772 


0.01620 


12.88 ±0.01 


0.82 ± 0.01 


SabA 


20 


10.87 ±0.19 


8.74 ± 0.03 


0.42 ± 0.01 


NGC5784 


0.01815 


12.85 ±0.01 


0.86 ± 0.01 


S0A 


1 


11.03 ±0.19 






UGC09598 


0.01857 


13.94 ± 0.02 


0.74 ± 0.03 


SbcAB 


14 


10.45 ±0.18 


8.67 ± 0.06 


0.27 ± 0.01 


UGC09665 


0.01838 


13.84 ± 0.02 


0.78 ± 0.02 


SbA 


7 


9.80 ±0.18 


8.63 ± 0.03 


0.03 ± 0.01 


NGC5829 


0.01865 


13.86 ± 0.02 


0.59 ± 0.03 


ScA 


24 


10.31 ±0.16 


8.53 ±0.10 


0.55 ± 0.02 


NGC5876 


0.01093 


12.87 ± 0.01 


0.91 ±0.01 


SOaB 





10.55 ± 0.20 






NGC5888 


0.02902 


13.36 ± 0.01 


0.86 ± 0.02 


SbB 





11.26 ±0.19 






UGC09777 


0.01547 


14.11 ±0.02 


0.72 ± 0.02 


SbcA 


8 


10.19 ±0.17 


8.65 ± 0.07 


0.33 ±0.01 


NGC5908 


0.01094 


12.48 ± 0.01 


1.02 ± 0.01 


SaA 


4 


10.83 ± 0.21 


8.67 ± 0.01 


0.53 ± 0.01 


NGC5930 


0.01861 


12.81 ±0.01 


0.85 ± 0.01 


SabAB 


5 


10.31 ±0.19 


8.68 ± 0.05 


0.32 ± 0.01 


UGC09873 


0.01844 


14.82 ± 0.04 


0.71 ±0.05 


SbA 


5 


10.05 ±0.17 


8.61 ± 0.06 


0.35 ± 0.01 


UGC09892 


0.01886 


14.46 ± 0.03 


0.72 ± 0.04 


SbcA 


7 


10.20 ±0.17 


8.67 ± 0.06 


0.25 ± 0.01 


NGC5953 


0.01671 


12.64 ± 0.01 


0.84 ± 0.01 


SaA 


4 


10.10 ±0.19 


8.70 ±0.12 


0.10 ±0.01 


ARP220 


0.01813 


13.31 ±0.01 


0.77 ± 0.02 


SdA 





10.74 ±0.18 






NGC5957 


0.01600 


12.72 ± 0.01 


0.74 ± 0.01 


SbB 


31 


9.95 ±0.18 


8.73 ± 0.08 


-0.27 ± 0.01 


IC4566 


0.01846 


13.45 ± 0.01 


0.87 ± 0.02 


SbB 


6 


10.80 ±0.19 


8.74 ± 0.07 


0.38 ± 0.01 


NGC5987 


0.01992 


12.35 ± 0.01 


0.96 ± 0.01 


SaA 





10.73 ± 0.20 






UGC 10043 


0.01715 


14.74 ± 0.07 


0.91 ± 0.09 


SabAB 


4 


9.45 ± 0.20 


8.58 ± 0.03 


-0.66 ± 0.01 


NGC6004 


0.01272 


12.91 ±0.01 


0.79 ± 0.02 


SbcB 


38 


10.56 ±0.18 


8.78 ± 0.03 


0.52 ± 0.01 


IC1151 


0.01712 


13.25 ± 0.01 


0.52 ± 0.01 


ScdB 


22 


9.62 ±0.15 


8.50 ± 0.06 


-0.02 ± 0.01 


UGC 10 123 


0.01239 


13.98 ± 0.02 


0.93 ± 0.02 


SabA 


3 


10.27 ± 0.20 


8.68 ± 0.03 


0.43 ± 0.01 


NGC6032 


0.01424 


13.58 ± 0.02 


0.83 ± 0.02 


SbcB 


2 


10.47 ±0.19 
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S. F. Sanchez et al.: The Mass-Metallicity relation explored with CALIFA: 
Table 1. Properties of the sample of galaxies. Continue 



Galaxy 


redshift 


V-band (mag) 


B-V 


Type 


Nhii 


log(Mass/MO) 


12+log(0/H) 


logCSFR/MOyr 1 ) 


NGC6060 


0.01446 


12.74 ± 0.01 


0.82 ± 0.01 


SbA 


16 


10.81 ±0.19 


8.69 ± 0.08 


0.86 ±0.01 


UGC 10205 


0.02185 


13.36 ±0.01 


0.88 ± 0.01 


SOaA 


3 


10.99 ±0.19 


8.63 ± 0.01 


0.52 ± 0.01 


NGC6063 


0.01938 


13.34 ±0.01 


0.65 ± 0.02 


SbcA 


38 


9.96 ±0.16 


8.61 ±0.10 


-0.02 ± 0.01 


IC 1199 


0.01557 


13.52 ±0.01 


0.73 ± 0.02 


SbAB 


17 


10.45 ±0.18 


8.77 ± 0.05 


0.42 ± 0.01 


NGC6081 


0.01679 


13.23 ±0.01 


0.94 ± 0.01 


SOaA 





10.87 ± 0.20 






UGC 10297 


0.01757 


14.34 ± 0.02 


0.60 ± 0.03 


ScA 


5 


9.26 ±0.16 


8.41 ± 0.09 


-0.35 ± 0.01 


UGC10331 


0.01589 


14.22 ± 0.03 


0.54 ± 0.03 


ScAB 


9 


9.92 ±0.15 


8.48 ± 0.08 


0.63 ± 0.01 


NGC6154 


0.01993 


13.40 ± 0.01 


0.82 ± 0.02 


SabB 


15 


10.81 ±0.19 


8.68 ± 0.04 


0.50 ± 0.01 


NGC6155 


0.01802 


12.87 ± 0.01 


0.64 ± 0.01 


ScA 


32 


10.03 ±0.16 


8.72 ± 0.05 


0.42 ± 0.01 


UGC10384 


0.01653 


14.25 ± 0.02 


0.75 ± 0.02 


SbA 


1 


10.21 ±0.18 






NGC6168 


0.01833 


14.27 ± 0.02 


0.65 ± 0.02 


ScAB 


11 


9.46 ±0.16 


8.51 ±0.05 


0.14 ±0.01 


NGC6186 


0.01957 


12.89 ± 0.01 


0.83 ±0.01 


SbB 


9 


10.36 ±0.19 


8.78 ± 0.04 


0.48 ± 0.01 


UGC 10650 


0.01097 


15.04 ± 0.03 


0.58 ± 0.03 


ScdA 


5 


9.27 ±0.16 


8.40 ± 0.02 


0.06 ± 0.01 


UGC10710 


0.02785 


14.05 ± 0.02 


0.82 ± 0.02 


SbA 





10.87 ±0.19 






NGC6310 


0.01128 


13.19 ±0.01 


0.82 ± 0.01 


SbA 


1 


10.39 ±0.19 






UGC 10796 


0.01040 


14.45 ± 0.03 


0.47 ± 0.03 


ScdAB 


4 


9.39 ±0.14 


8.44 ± 0.07 


-0.30 ± 0.02 


NGC6361 


0.01267 


13.06 ± 0.02 


0.99 ± 0.02 


SabA 


6 


10.73 ± 0.21 


8.68 ± 0.04 


0.98 ± 0.01 


UGC10811 


0.02905 


14.16 ±0.02 


0.82 ± 0.03 


SbB 





10.89 ±0.19 






IC1256 


0.01568 


13.60 ± 0.01 


0.72 ± 0.01 


SbAB 


25 


10.40 ±0.17 


8.76 ± 0.05 


0.42 ± 0.01 


NGC6394 


0.02879 


14.14 ±0.02 


0.85 ± 0.03 


SbcB 





10.89 ±0.19 






UGC 10905 


0.02580 


13.31 ±0.01 


0.95 ± 0.02 


SOaA 





11.25 ±0.20 






NGC6478 


0.02234 


13.76 ±0.01 


0.82 ± 0.02 


ScA 


30 


10.79 ±0.19 


8.67 ± 0.05 


1.06 ±0.01 


NGC6497 


0.02013 


13.26 ±0.01 


0.82 ± 0.02 


SabB 


9 


10.91 ±0.19 


8.70 ± 0.01 


0.61 ±0.01 


UGC 11262 


0.01836 


14.81 ± 0.04 


0.59 ± 0.05 


ScA 


8 


9.91 ±0.16 


8.62 ± 0.04 


-0.05 ± 0.01 


NGC6978 


0.01986 


13.15 ±0.02 


0.84 ± 0.02 


SbAB 


9 


10.94 ±0.19 


8.73 ± 0.04 


0.60 ± 0.01 


UGC 11 649 


0.01249 


13.22 ± 0.03 


0.85 ± 0.03 


SabB 


18 


10.50 ±0.19 


8.69 ± 0.09 


-0.13 ±0.02 


UGC11680NED01 


0.02649 


13.86 ± 0.02 


0.90 ± 0.03 


SbB 


10 


11.01 ±0.20 


8.68 ± 0.06 


0.87 ± 0.01 


NGC7047 


0.01907 


13.41 ± 0.02 


0.83 ± 0.02 


SbcB 





10.82 ±0.19 






UGC11717 


0.02128 


13.90 ± 0.02 


1.03 ± 0.03 


SabA 


2 


10.94 ± 0.21 






MCG-01-54-016 


0.01062 


14.83 ± 0.05 


0.42 ± 0.06 


ScdA 


1 


9.17 ±0.14 






UGC 11740 


0.02118 


14.26 ± 0.03 


0.71 ±0.04 


SbcA 


9 


10.42 ±0.17 


8.65 ± 0.03 


0.57 ±0.01 


UGC 11792 


0.01567 


14.49 ± 0.04 


0.75 ± 0.05 


SbcA 


4 


10.06 ±0.18 


8.68 ± 0.01 


0.44 ± 0.01 


UGC 12054 


0.01674 


14.42 ± 0.03 


0.51 ±0.04 


ScA 


5 


8.97 ±0.15 


8.32 ± 0.05 


-0.59 ± 0.01 


NGC7311 


0.01503 


12.19 ±0.01 


0.88 ± 0.01 


SaA 


10 


11.17 ±0.19 


8.75 ± 0.02 


0.72 ± 0.01 


NGC7321 


0.02362 


13.16 ±0.01 


0.75 ± 0.02 


SbcB 


37 


11.01 ±0.18 


8.71 ± 0.05 


0.98 ± 0.01 


UGC12185 


0.02184 


13.82 ± 0.02 


0.76 ± 0.02 


SbB 


6 


10.66 ±0.18 


8.66 ± 0.01 


0.29 ± 0.01 


UGC 12224 


0.01171 


13.66 ± 0.03 


0.84 ± 0.04 


ScA 


48 


10.23 ±0.19 


8.64 ± 0.08 


0.07 ± 0.01 


UGC 12308 


0.01724 


14.56 ± 0.04 


0.45 ± 0.05 


ScdA 


10 


8.99 ±0.14 


8.30 ± 0.08 


-0.52 ± 0.01 


NGC7466 


0.02489 


13.80 ± 0.02 


0.75 ± 0.02 


SbcA 


20 


10.77 ±0.18 


8.66 ± 0.03 


0.84 ± 0.01 


NGC7489 


0.02072 


13.37 ± 0.02 


0.55 ± 0.02 


SbcA 


24 


10.56 ±0.15 


8.52 ±0.13 


1.08 ±0.01 


NGC7536 


0.01540 


13.42 ± 0.02 


0.67 ± 0.03 


ScAB 


29 


10.42 ±0.17 


8.61 ± 0.08 


0.70 ± 0.01 


NGC7549 


0.01549 
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Appendix A: Dependence of the ai-z relation with 
the SFR: Apperture effects. 

We show in the main text our results on the study of the pos- 
sible dependence of the AI-Z relation with the SFR. Using the 
IFS data provided by the CALIFA survey, we were able to de- 
rive the AI-Z relation with abundances not affected by aperture 
effects. In a similar way we derive the complete integrated SFR 
for t hese galaxies. Contrary to previous results in t he literature 
(e.g. iMannucc i et alj 120101 : iLara-Lopez et al.ll2010l) we find no 
secondary dependence of the abundance with the SFR once con- 
sidered the dependence of both parameters with the Mass. The 
only similar result found in the literature where no dependence 
on the SFR is described used a sam ple of integrated spe ctra taken 
using drift-scan observations (e.g. Hughes et al. 2012). In sum- 
mary, when the biases due to aperture are minimized the depen- 
dence on the SFR seems to disappear. 

It is beyond the scope of this article to explain why other 
studies find a correlation that we cannot reproduce, in particu- 
lar when we use a less biased and more detailed estimation of 
the considered quantities. However, as a sanity check, and in or- 
der to provide a suitable explanation, we performed a test of the 
aperture effects. The oxygen abundance and the SFR are poten- 
tially affected by this aperture effect. However, the stellar mass is 



the less affected, in particular if only photometric data is consid- 
ered. With CALIFA we are in the unique situation of being able 
to reproduce the aperture effects present in single aperture obser- 
vations, by simulating the considered apertures and re-analysing 
the corresponding simulated spectra. 

It is important to use as our starting dataset for the simu- 
lations a sample that we know a priori that it does not show 
any secondary correlation between the abundance and the SFR, 
when the dependency with the Mass is removed. In this case we 
will be sure that the induced correlation was not present as an 
input ingredient of the simulation. 

For each galaxy in our dataset with well defined character- 
istic oxygen abundance, we extracted a set of aperture spectra 
centred in the peak emission of the galaxy and with consecu- 
tive apertures corresponding to the distance of each of the de- 
tected H ii regions from our catalogue. This ensures that there 
is a clear change in the ionization conditions between consec- 
utive apertures. Each of these aperture spectra correspond to 
the same galaxy observed at different r edshifts. To be consi stent 
with the redshift ranges adopted by Mannucci et al. (2010) and 
lLara-Lopez et a l. (2010), we restricted our apertures to diame- 
ters larger than 12". This roughly corresponds to the projected 
size of a SDSS fiber at the lower redshift adopted by this study 
(z ~0.07). Once extracted the spectra, we repeat the same anal- 
ysis that we perform for the individual H n regions, described 
in Sec. 13. II deriving the oxygen abundance and integrated SFR. 
Together with the integrated stellar masses, we ended up with 
~2500 individual simulations of the three quantities analysed in 
this study. 

Finally, we perform a Monte-Carlo simulation to take into 
account the possible effects of the noise in the overall process. 
For each of the three considered quantities 20 individual realiza- 
tions were derived, taking into account the corresponding propa- 
gated errors (including also those of the stellar masses). We end 
up with ~50,000 simulated aperture biased quantities. 

Figure IA.1I shows the results of this simulation. In the top- 
left panel we show the derived AI-Z relation for the simulated 
spectra. As expected it follows the observed relation, presented 
in Fig. [4] This result was already shown by Rosales-Ortega etaD 
(2012), where they reproduce the global A4-Z from the local one 
(derived from the abundance gradients), based purely on simu- 
lated galaxies. However, a more detailed comparison shows that 
although both relations are very similar there are slightly dif- 
ferences: (i) the asymptotic oxygen abundance for large masses 
is ~0.05 dex larger in the simulated relation and (ii) the dis- 
persion of oxygen abundances along the average value is larger 
(o"iog(o/H) ~0. 1 dex, similar to the one reported by T04). Both re- 
sults are expected: (1) we are including more central estimations 
of the abundance for the same range of masses. In the center 
of the galaxies th e abundance is larger than in the outer parts 
by ~0.1 dex (e.g. Sanc hez et alJl2012bl) . Thus, the net effect is 
an increase of the upper envelop of the abundance; (2) we are 
including measurements of different abundances for the same 
mass range as we enlarge the covered area of the galaxy, sim- 
ulating the effect of the redshift. Therefore, the dispersion at a 
certain mass range is increased. 

Figure [ATI top-righ panel shows the distribution of the sim- 
ulated SFRs along the stellar masses. As expected we reproduce 
the relation between both quantities, in a similar way as de- 
scribed in Fig. |4] The main differences are (1) the dispersion is 
slightly larger: cr log(SFR) -0.35 dex, compared to cr log(SFR) ~0.18 
dex; (2) there is a tail to lower values of the SFR at a given mass. 
This is in particular evident for the low-mass range of the figure 
(M ~ 10 95 Mo). Once more, these effects are clearly explained 
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Fig. A.l. Top-Right panel: Distribution of the oxygen abundance along the stellar mass for the -50,000 simulated aperture mea- 
surements, as described in the text. The contours indicate the same encircled fraction of object than in precendent figures (e.g., 
Fig. 0. The solid yellow circles indicate the average abundances for consecutive bins of 0.2 dex in stellar mass. Only those binned 
values with a significant number of objects are shown. Top-left panel: Distribution of the integrated star-formation rate along the 
stellar mass for the simulated galaxies. Contours and symbols represent the same as in previous panel. The solid line shows the 
regression line between the two parameters. Bottom-left panel: Distribution of the differential oxygen abundance with respect to 
the derived M-Z relation along the SFR for the simulated galaxies, when only the aperture effect is considered. Bottom- right panel: 
Distribution of the differential oxygen abundance with respect to the derived M-Z relation along the SFR for the simulated galaxies, 
when both the aperture effect and the evolution of the M-Z relation with redshift is taken into account. In both panels the contours 
and sybmols represent similar concepts like in previous figures and panels. The solid line shows the best derived fit using a simple 
linear regression method. 



by the aperture bias, since at a given mass we are sampling ar- 
eas of lower integrated SFR than the real one (i.e., the central 
regions), which is an upper envelop. 

Finally, we derive the differential oxygen abundance with re- 
spect to the considered M-Z relation, for the simulated data. We 
try to determine if there is any correlation with the SFR, fol- 
lowing the analysis in Sec. 13.31 Figure IA.1I bottom-left panel 
shows the distribution of these two parameters. There is a clear 
correlation (r - 0-92 ), in the simil ar sense as d escribed by 
iMannucci et alJ d2010h andlLara-Lope zetal] d2010l) . i.e., galax- 
ies with lower SFR seem to be more metal poor. Once applied 
this correlation, we recover the initial dispersion of the M-Z rela- 
tion described for our original dataset (cr Alog(0/H) ~0.07 dex), 
a decrease similar to the one reported by these two authors. The 
dependence of the SFR is Ai og (o/H) °c— 0. 1 log(SF R), just slightly 
weaker than the linear term of the equation 2 of IMannucci et al.l 



(2010), -0.14 dex/log(SFR), or the combination of equations 3 
and 4 of the same authors. 

In order to investigate further this issue, we have included 
in our simulati on the reported M-Z evolution with redshift, as 
described by Moustakas et al.l (1201 ll) . For doing so, we have de- 
termined at which redshift corresponds each of the simulated 
apertures, and we apply the described correction to the corre- 
sponding derived abundance. Results are shown in Figure IA.1I 
bottom-right panel. Taking into account the possible redshift de- 
pendence, the new differential oxygen abundance have a sharper 
dependence on the SFR, Alog(O/H)oc-0.21og(SFR). 

Despite of the fact that this is a very simplistic simulation, 
based on a reduced number of input galaxies, in general we 
have shown that aperture effects, together with other evolution- 
ary properties of galaxies, may induce the secondary correlations 
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between the A1-Z and the SFR that we are not able to reproduce 
using our IFS data. 
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